Possible Solution of the J/i/' Production Puzzle 
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We argue that the s-channel cut contribution to J/i/; hadroproduction can be significantly larger 
than the usual cut contribution of the color-singlet mechanism (CSM), which is known to under- 
estimate the experimental measurements. A scenario accounting for intermediate cc interactions is 
proposed that reproduces the data at low- and mid-range transverse momenta Pt from the Fermilab 
Tevatron and BNL Relativistic Heavy Ion Collider. The J/i/; produced in this manner are polarized 
predominantly longitudinally. 

PACS numbers: 13.60.Le, 11.40.-q, 13.85.Ni, 14.40.Gx PRL 100, 032006 (2008)— arXiv:0709.3471v2 



Although heavy quarkonia are among the most ana- 
lyzed bound quark systems, ever since the first measure- 
ments by the CDF Collaboration of the direct production 
of J/ip and i})' at ^/s — 1.8 TeV we are facing persis- 
tent disagreements between theoretical predictions from 
various models and experimental studies of the cross sec- 
tion and the polarization (for reviews see [§]). 

The recent confirmation by CDF [3] of their previous 
polarization measurement 0] showing an unpolarized or 
slightly longitudinally polarized J/^ yield has reinforced 
doubts about the applicability of the quark-velocity ex- 
pansion [v) of nonrelativistic QCD (NRQCD) for the 
rather "light" cc system. On the theory side, many new 
results completed — but also questioned — our knowl- 
edge of charmonium production: the long-awaited next- 
to-leading order (NLO) QCD corrections to the color- 
singlet contributions are now available and show sig- 
nificant enhancement of the cross section; an up-to- 
date proof [3] of NRQCD factorization holding true at 
any order in v in the gluon-fragmentation channels was 
provided; the universality of the nonpcrturbative in- 
put of NRQCD was challenged by fixed-target measure- 
ments ^ , similar to what had been found previously for 
the HERA data 0]; NRQCD factorization was shown to 
require modification in fragmentation regions where three 
heavy quarks have similar momenta [lO|; finally, the c- 
quark fragmentation approximation was shown [ll| to 
be only valid at much higher Pt than expected from the 
pioneering works of Ref . . 

In view of the difficulty of most theoretical approaches 
in re pro ducing the experimental data, the authors of 
Ref. 3| undertook a systematic study of the cut con- 
tributions due to off-shell and nonstatic quarks. In par- 
ticular, they questioned the assumption of the CSM that 
takes the heav y q uarks forming the quarkonium (Q) as 
being on-shell If they are not, the usual s-channel 
cut contributes to the imaginary part of the amplitude 
and need to be considered on the same footing as the 
CSM cut. These s-channel cut contributions, which were 
the specific focus in Ref. , are illustrated in Fig. [1] 

In order to provide a conserved current for such off- 



shell configurations, one must introduce an additional 
four-point function, or contact current [l5| . Dynamically 
such a current accounts for the interactions between the 
cc pair emitting the external gluon. In fact, this mech- 
anism arises because of the possibility that the outgoing 
gluon stems from within the dressed cc J/ ip vertex, as de- 
picted in Fig. [21 The quark pair (ci, C2) that makes the 
final J / i}}-g\non state is now in a color-octet state which 
thus recovers the necessity for such configurations as a 
natural consequence of restoring gauge invariance. 

The construction of this 4-point function based on 
gauge invariance alone is not unique since additional 
transverse contributions have no bearing on current con- 
servation. In principle, the details of this 4-point func- 
tion and of the necessary transverse contributions would 
follow from a full dynamical treatment that consistently 
accounts for all such interactions. At present at least, 
this is beyond the scope of what is theoretically possible. 

Yet, it is possible to construct a 4-point function that 
satisfies certain minimal requirements U 

Ei. The 4- 

point function proposed in ^] provided a conserved cur- 
rent but was not entirely satisfactory since it contained 
poles (similar to the basic direct and crossed contribu- 
tions), and such poles for the contact current are un- 
physical and therefore should be avoided [f^. 

To this end, drawing on the experien ce g ained in pion 
photoproduction processes E H Eof, we propose 
here a contact current for J / tp production which satisfies 
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FIG. 1: (Color online). (a,b) Leading-order (LO) s-channel 
cut diagrams contributing to gg — > Qg with direct and crossed 
box diagrams employing the cc J/tp vertex. The crosses indi- 
cate that the quarks are on-shell. (c) Box diagram with ccQg 
contact term mandated by gauge invariance. 
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FIG. 2: Example of the mechanisms (a) contributing to the 
dressing of the 3-point function r'"^-* and (b) responsible for 
the 4-point function F'*': the external gluon is attached here 
within gluon loops of the dressed vertex, thus producing a di- 
agram without poles and with a kinematic behavior genuinely 
different from the initial 3-point vertex. 



the requirements of full gauge invariance, going beyond 
mere current conservation, in terms of the generalized 
Ward-Takahashi identity [3, [l3| ■ The basic mechanism 
how to do this was proposed a long time ago by Drell 
and Lee [15| for the electromagnetic case, where it basi- 
cally amounts to the minimal substitution prescription 
^/^ ^ + iQA^' (Q: charge, A^'■. vector potential) 
in an effective Lagrangian corresponding to the dressed 
hadronic vertex. As noted by Drell and Lee, this par- 
ticular approach is deficient in that it violates the high- 
energy scaling behavior, because in avoiding poles for 
the 4-point function, it partially replaces the true mo- 
mentum dependence of the vertices by constants. We 
avoid this shortcoming by building a phenomenological 
vertex incorporating this low-energy behavior and the 
expected momentum dependence of the vertices at high 
energies by interpolating between these two regimes. As 
we show below, the resulting enhanced s-channel cut con- 
tributions reproduce the experimental data from both the 
CDF and PHENIX collaborations up to transverse mo- 
menta Pt ^ 10 GeV. 

The present formalism builds on the approach devel- 
oped in [l3| . where the transition qq —> Q is described 
by the 3-point function 



r(f)(p,p) = r(p,p)7^ 



(1) 



where P = pi — p2 and p = {pi + p-i)l^ are the total 
and relative momenta, respectively, of the two quarks 
bound as a quarkonium state, with p\ and pi being their 
individual four-momenta. Ansatz ([T]) amounts to rep- 
resenting the vector meson as a massive photon with 
a nonlocal coupling. In the present work, we describe 
the relative- momentum distribution r(p, P) of the quarks 
phenomenologically as a Gaussian, as in Ref. where 
full details can be found. 

The generic picture of the physical origin of the dressed 
vertex function r(p, P) is given in Fig. [2{a), and the en- 
suing contact current when coupling the gauge boson into 
this vertex is illustrated in Fig.[2ljb). 

The requirement of gauge invariance is usually written 
in terms of the generalized Ward-Takahashi relations [l^ 
for the complete current. For the purpose of restoring 



gauge invariance, it is most convenient to rewrite this into 
an equivalent condition for the contact current (see [l7j 
for the analogous procedure in pion photoproduction) . 
To this end, let us write the 4-point function depicted in 
Fig. mb) as 



r(4) = -^g,Tt^M:^^ 



(2) 



where gs is the strong coupling constant, T-"^ the color 
matrix, and /i and v are the Lorentz indices of the outgo- 
ing J/i/i and gluon, respectively. For simplicity, we have 
suppressed all indices on the left-hand side. The ccJ/tjj 
vertex function F^^-* with the kinematics of the direct 
graph is denoted here by Fi and for the crossed graph by 
Fs, I.e., Fi =F(ci-f,P) and F2 = F (c2 4- f , P) , as 
shown in Figs.[IJa) and[IJb). The gauge-invariance con- 
dition for the contact current for the outgoing gluon 
with momentum q reads now 



q,Mc = Ti - F2 



(3) 



since this is precisely the four-divergence contribution 
needed to cancel the corresponding terms arising from 
the four-divergences of Figs. [Ija) and[ljb). 

We emphasize here that, within the present semi- 
phcnomcnological approach, the procedure to preserve 
QCD gauge invariance for the gluon coupling follows ex- 
actly along the lines of QED for a photon since the four- 
point function ([2]) factorizes in the color matrix and the 
gluon coupling. This finding is directly related to the fact 
that the two vertex functions on the right-hand side of 
the four-divergence contribution ([3]) arise from the gluon 
coupling to the two intermediate quark lines in Figs.[TJa) 
and[T]Jb) which, apart from color factors, is exactly like 
a photon coupling to spin-1/2 particles. 

We now employ the usual construction [H, [l^l 
for the contact current in terms of an auxiliary function 
F — P(ci, C2, q) and put 



Mr 



{2c2 + qY [Ti-F) , (2ci - g)'' (F2 - P) 



(C2 + 



m, 



(ci — q)'^ — m? 



(4) 



where we take c\ 



m? and P^ = A'P from the 



beginning, with m and M being the masses of the quark 
and the J/'i/'i respectively. One easily verifies that this 
additional contact current satisfies the gauge-invariance 
condition ([3]). It is found, in particular, that P cancels 
out in the four-divergence. 

The function P(ci,C2,(i') must be chosen so that the 
current ([4]) satisfies crossing symmetry (i.e., symmetry 
under the exchange ci <-> — C2) and is free of singulari- 
ties. The latter constraint implies P = Fq at either pole 
position, i.e., when (c2 -I- 9)^ — rri^ or (ci — g)^ = m^, 
where the constant Fq is the (unphysical) value of the 
momentum distribution F(p, P) when all three legs of the 
vertex are on their respective mass shells. In principle, 
employing gauge invariance as the only constraint, we 
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may take F = Tq everywhere. This corresponds to the 
minimal substitution discussed by Drell and Lee [13] (for 
a complete derivation see [2l|) who pointed out, however, 
that this does not provide the correct scaling properties 
at large energies, which means within the present con- 
text that F = Tq would not lead to the expected Pt 
scaling of the amplitude. Numerically, this choice over- 
shoots the experimental data by more than one order of 
magnitude at Pt = 20 GcV. In fact, in the large relative- 
momentum region, we expect the contact term and there- 
fore the function F{ci, C2, q) to exhibit a fall-off similar to 
the vertex functions themselves, contrary to the minimal 
substitution procedure. The simplest crossing-symmetric 
choice with this behavior is F = Ti -I- - riTz/ro [H . 
The solution we propose here is to build F(ci, C2, q) from 
these two limiting cases. To this end, it is natural to 
choose the following simple ansatz 



F{ci,C2,q) = To 



, (ro-ri)(ro-r2) 

- n\ciC2} , 

i 



(5) 



where the (crossing-symmetric) function h{ciC2) rises to 
become unity for large relative momentum. Note that it 
is not necessary that h{ciC2) actually vanish at either pole 
since the factor on its right vanishes for (c2 -I- q)^ — m? 
or (ci — g)^ = m? and we recover the choice of = To at 
either pole. The phenomenological choice for the inter- 
polating function h{ciC2) used in our calculations is 



h{ciC2) = 1 — a 



— (ciC2 + TO^) 



(6) 



with two parameters, a and n. We would like to empha- 
size at this point that this choice is in no way unique. In a 
manner of speaking, this choice is simply a way of param- 
eterizing our ignorance by employing minimal properties 
of r^^^. We shall not, however, discuss other choices here 
since our main motivation is to show that s-channel cut 
contributions can be large and can indeed reproduce the 
data — the physical picture that emerges could then be 
tested in other production regimes. 

In the kinematical region accessed at the Tevatron, the 
direct J/ 1/' are produced by gluon fusion and a final-state 
gluon emission is required to conserve C-parity and pro- 
vide the J jij] with its Pt- The relevant diagrams for the 
LO gluon fusion process can be found in [l3|, the only 
difference in the present treatment being the new choice 
ofrW. Also, we use the same normalization of as 

in 

The double-differential polarized cross section in trans- 
verse momentum Pt and rapidity y is given by [l3| 



dydPT 



dxi 



2sPTg{xi)g {x2{xi)) dor 
Vi(yia:i - ETeV) ~d 



(7) 



where dar/dt is the partonic differential cross section, 
with r = L,Ti,T2 being the quarkonium helicity, and 



s = (fci -|- ^2)^, i = {k2 — qY and it = {ki — q)"^ are the 
Mandelstam variables for the partonic process. In the 
present calculations, we use the LO gluon distribution 
g{x) of 22], and the same mass and size parameter A as 
given in Ref. [l^. In any case, our conclusion that the s- 
channel cut contribution can reproduce the experimental 
data with adjustments of the values for a and k would 
not change at all over a wide range for these parameters. 

Figure [3Ka) shows our results with parameter val- 
ues a = 4 and K = 4.5 GeV for ^ = 1.8 TeV in 
the pseudorapidity range \ri\ < 0.6 compared with the 
cross-section measurement of direct J/ip by CDF 0], 
the usual LO GSM from go ^ J/V' 5 Q and LO 



GSM from gg — > J/tJjcc 11|. Our results agree very 
well with the GDF data up to about Pt = 10 GeV. 
At higher Pt, our curve falls below the data as ex- 
pected from the genuine 1/Pt scaling of a LO box dia- 
gram. Inclusion of higher-order corrections incorporating 
fragmentation- type topologies (~ ^/Pt) and associated- 
production channels are expected to fill the gap between 
data and theory at high Pt- It is interesting to note 
the different Pt behaviors of ctt and ctl leading to a 
dominance of the latter at large Pt and a negative value 
for the polarization a. Figure [SJb) shows our results at 
= 200 GeV, still with a = 4 and k = 4.5 GeV, com- 



LO CSM: J/v+cc 
LO CSM: J/y+g 
CDF data 




FIG. 3: (Color online), (a) Comparison between polarized 
((tt and ctl) and unpolarized (crtot) cross sections [with pa- 
rameters a = 4, K = 4.5 GeV in Eq. ®], LO CSM con- 
tributions, and CDF experimental data [3] at the Tevatron 
(-yi = 1.8 TeV, pseudorapidity \ri\ < 0.6). (b) Compari- 
son between ctt, o"l, o"tot and PHENIX data [23] at RHIC 
{y/s = 200 GeV, rapidity \y\ < 0.35). 
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FIG. 4: Prompt J/ip polarization: theory vs. CDF data 



pared with the PHENIX data [23 1. 

Since J/?/; polarization measurements exist only for the 
prompt yield, we have computed a from our direct- J/i/' 
cross sections in two extreme cases, one where the J/ip^s 
from Xc are 100% transverse and another where they are 
100% longitudinal, the first scenario being the more likely 
one. Figure |4] shows the comparison between this com- 
putation and the recent results by CDF at ^/s — 1.96 
TeV A]. 

In conclusion, we have shown here that among the two 
singularities in the box diagram contribution to quarko- 
nium production at LO, the one from the s-channel cut 
is significantly larger than expected when one includes 
interactions between the quark pair which binds into 
the quarkonium and emits the final-state gluon. In the 
present work, we have accounted for the effect of those 
interactions by a phenomenological contact current, or 
4-point function, F^^^ which follows directly from the im- 
plementation of the full gauge-invariance requirements 
appropriate for a dressed vertex F^^^ . This 4-point func- 
tion, however, is not fully constrained by gauge invari- 
ance which permitted us to interpolate between the min- 
imal substitution discussed earlier by Drell and Lee and 
the expected large relative-momentum behavior. 

As we showed, this 4-point function provides a much 
larger contribution than the direct and crossed ones con- 
taining 3-point functions (cf. Fig.[T]), and this can easily 
bring about agreement with the experimental data. In 
NRQCD, color-octet matrix elements account for transi- 
tions between a colored heavy-quark pair into a quarko- 
nium by soft unseen gluon emission in the final state. 
In the present approach, the 4-point function accounts 
for gluon exchanges between the heavy quarks which 
emit the final-state gluon. As for the matrix elements 
of NRQCD, which are unknown and then fit, we fixed 
the unconstrained parameter of this function in order to 
reproduce the experimental data at ^/s = 1.8 TeV from 
the CDF collaboration at the Tevatron for Pt < lOGeV. 

With the same parameters, we also obtain a very 
good description of the experimental measurements from 
PHENIX at RHIC at ^/s ^ 200 GeV. Moreover, our 
prediction for the polarization for the prompt J/ip yield 
is mostly longitudinal. This looks promising since one 
expects contributions from the CSM cuts — which are 
known to be enhanced by the NLO corrections 0] — and 



from the real part which is not evaluated at present time. 

A similar enhancement by inclusion of the s-channel 
cut is expected in all production processes where the 
J/^p is associated with a gluon, e.g., photon-photon col- 
lision at LEP as well as in photo- and lepto-production 
at HERA. 
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